Stability of Solid State Reaction Fronts 
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We analyze the stability of a planar solid-solid interface at which a chemical reaction occurs. 
Examples include oxidation, nitridation, or silicide formation. Using a continuum model, including a 
general formula for the stress-dependence of the reaction rate, we show that stress effects can render 
a planar interface dynamically unstable with respect to perturbations of intermediate wavelength. 



Many important reactions occur at solid-solid inter- 
faces, and require one of the reacting species to diffuse 
to the interface through one of the solids. The oxidation 
of silicon, a complex process vital to the fabrication and 
function of silicon devices, is the best known example; 
but the oxidation of metals, nitridation of silicon, and sili- 
cide formation fall into the same general category. Since 
the original and reacted materials typically have different 
lattice constants, the reaction generates stress, which in 
turn alters the reaction rate. In non-planar geometries, 
where the stress is non-uniform, the resulting structure 
can be drastically affected . There is a longstanding 
effort to understand these stress effects on morphology 
in Si oxidation and in other systems 

It is well known that in the case of epitaxial growth 
of a solid at a free surface, stress causes a morphological 
instability in an initially planar surface j^j]. In this paper, 
we show that the stress at a solid-solid reaction front can 
similarly lead to an instability Q . 

The instability here, however, differs in two impor- 
tant respects from that occurring in epitaxial growth. 
First, the instability in epitaxial growth is driven by the 
thermodynamics, always acting to lower the energy of 
the system. In contrast, the reaction instability is es- 
sentially dynamical in nature. It results from the effect 
of stress on the reaction rate rather than on energetics. 
Stress may either stabilize or destabilize the planar re- 
action front, even though the energy is always lowered 
by long-wavelength deviations from planarity. Second, a 
free surface under stress is always unstable at long wave- 
lengths and stable at short wavelengths. In contrast, 
solid-solid reaction fronts are stabilized at long and short 
wavelengths by diffusion and interface-tension effects, re- 
spectively. An instability can therefore occur only at 
intermediate wavelengths. 

The system we study is illustrated schematically in Fig. 
1. For simplicity, we use the language of oxidation to de- 
scribe it. However, the model we now discuss is quite 
general, and applies equally well to the other reactions 
mentioned above. In oxidation, the surface layer of a solid 
(typically metal or semiconductor) is in contact with a 
reservoir of oxidant, such as O2 gas. The surface oxi- 
dizes, producing a thin solid film of oxide, as shown in 
Fig. 1. The oxidant must diffuse through the film to the 
film-substrate interface in order for further oxidation and 
film growth to occur. 



Typically there is some volume change upon oxidation, 
which produces stress 0. This in turn affects the oxi- 
dation and diffusion rates. In some materials, such as 
Si02, viscous flow of the oxide film or other inelastic pro- 
cesses can relieve these stresses, in whole or in part. Deal 
and Grove Q proposed a simple model for the oxida- 
tion of planar substrates in the high-temperature regime, 
where viscous flow is so rapid that stress relaxation can 
be treated as instantaneous. Their model has become 
the standard framework within which many subsequent 
oxidation problems involving curved geometries Q have 
been analyzed. 

Here we address the opposite regime, where the system 
is purely elastic and there is no viscous flow. This case 
is important for oxidation and other reactions at lower 
temperatures. As device structures shrink toward the 
nanometer scale, and processing temperatures are cor- 
respondingly reduced, the stress-induced instability we 
describe may well become important in systems where it 
was previously irrelevant. 

Consider the geometry of Fig. 1, where the substrate- 
film interface is a plane modulated by a small sinu- 
soidal perturbation with amplitude a q and wavelength 
A = 2ir/q. Our task is to evaluate the interface veloc- 
ity, and determine whether the deviation from planarity 
grows with time. If for some q, a q grows as the reaction 
proceeds, then the planar interface is unstable. 

The velocity v of a point xi on the reaction front is 



v(xj) = R(xi)p f 1 n I = [c(xi)ki - Pcfo] P f 1 n I . (1) 



Here R is the number of oxidant molecules per unit area 
that react at the interface per unit time to form oxide, 
fii is the local unit normal to the interface (Fig. 1), pf is 
the density of bound oxidant molecules in the oxide, c(x) 
is the concentration of diffusing oxidant molecules in the 
film, k\ and hi are the respective rate constants for the 
forward and reverse reactions (oxidation and reduction) 
at the interface, and p c is the density of sites through 
which oxidant can diffuse in the film. The Fourier com- 
ponent of this velocity with wave number q gives directly 
the rate of change of a q , and hence determines the sta- 
bility of a planar interface. 
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FIG. 1. An illustration of the geometry of the solid state 
reaction. The diffusion of the mobile reactant is represented 
by a random walk. 



With the current 
molecules given by j 



density j of diffusing oxidant 
= — £>Vc , , c(x, t) satisfies the 



diffusion equation dc/dt = D\7 2 c, where D is the oxi- 
dant diffusion constant. Typically, the diffusion is rapid 
enough to maintain a quasi-steady-state concentration 
pOf , so c obeys simply 



V^c = 



(2) 



The boundary condition on this equation at the in- 
terface comes from the requirement that all oxidant 
molecules flowing into the interface react with the sub- 
strate to produce oxide. This can be written 



n i ' 3 = R{xi) = c(a;l)fci - p c k 2 



(3) 



The boundary condition at the upper surface is that the 
normal oxidant current equals the rate at which oxidant 
is incorporated into the upper surface of the film from 



the reservoir 10 



j(xs) = -h [c* - c(xs)] n s 



(4) 



Here ft. is a rate constant (Henry's constant), x*s is a point 
on the upper surface, c* [> c(xs)] is the concentration 
of oxidant in the oxide film in equilibrium with the reser- 
voir, and hs is the unit normal to the upper surface. 

The solution of Eq. (||) subject to the boundary con- 
ditions (||) and ([|) is complicated by the fact that the 
rate constants fci and k 2 depend on the local stress. The 
stress dependence has been discussed from a phenomeno- 
logical perspective by several authors in the context of 
silicon oxidation Jl| . Here we apply a recent more com- 
plete treatment Wu, which we now briefly sketch. The 
reaction proceeds from an initial state of free energy F\ 
to a final state of free energy F 2 , through a transition 
state (the saddle point of the energy surface) with free 
energy F t . The local reaction rates can then be written 
as: 



where = 1/ksT, ks is Boltzmann's constant, i = 1 
or 2, and fco is a rate constant reflecting the microscopic 
"attempt frequency." Fx, Ft, and F 2 represent free en- 
ergies coarse-grained over distances large compared to 
atomic dimensions but small compared to other dimen- 
sions of the system. In general, Fx, F t , and F 2 depend 
on the stresses and interface curvature throughout the 
system. Note that we are concerned with free-energy 
changes Fi — F t associated with a reaction event, and 
these depend only on the local stress and curvature. 

We are considering stability with respect to an in- 
finitesmal perturbation from planarity, so we expand the 
stress and reaction rate to lowest order in a q . Then 
Fi = Fi + 5Fi, for i = l,t,2. Here F, is the (local) 
free energy for the planar interface, including stress ef- 
fects, and 5Fi is the change due to the curvature and 
to the extra stress produced by the sinusoidal corruga- 
tion. Since SF is 0(a q /X), we can take SF <C feT, 
so h w fcj [1 — j3{SF t — SFi)], with i = 1,2 and fc, = 
/coe^^' - ^- 1 . The reaction rate is then 

R[Si) « R{xi) + /3c{xi)kx(SF t - SFi) 

-l3 Pc k 2 (5F t -SF 2 ) , (6) 

where R(xj) = \c(Si)kx — p c k 2 ~\ is the rate for the flat 
interface. 

Expanding SF to lowest order in the local strain and 
curvature gives 



SF t - SFi = (7i« + o-fSe nn + o-f5e pp )p f 1 



(7) 



Here ji is the curvature-derivative of F t — Fi, and acts 
as an effective interface tension; k is the local curvature 
of the interface; Se nn and 5e pp are the extra strains nor- 
mal and parallel to the interface due to deviations from 
planarity; and erf and af are coefficients reflecting the 
stress of the transition state relative to that of the ini- 
tial or final state in the normal and parallel directions 
respectively Jl^] . Note that the expansion for SFt — SFi 
can be written in terms of the strains Se nn and Se pp of 
either the film or the substrate; but the coefficients af 
and erf assume different values in the two cases. 

We evaluate the strains for the geometry of Fig. 1 using 
linear elastic theory Jl3[ . The volume expansion accom- 
panying oxidation gives rise to elastic displacements u(x) 
of the material at position x. Taking both the substrate 
and oxide to be isotropic elastic media but with differ- 
ent elastic constants, one expresses the stress tensor a in 
terms of u through the standard relation 



A «,/(V . u s -f)S l3 + p sJ (du^/dxj + duf/dxi) 



"$K sJ r) s} S 



(8) 



h = k e-^- F '- 



(5) 



Here rj is the misfit strain parallel to the surface, with 
rf = and — r\ for the substrate (s) and film 
(f), respectively; K s f are the bulk moduli, K s f = 



2 



A s,/ + 2fi s ^ /3; and (A s ^, ^ s ^) are the Lame coefficients. 
The above form of the stress tensor incorporates the re- 
quirement that the stress-free state for the substrate has 
zero displacement, u s — 0, while the stress-free state for 
the film is achieved through a uniform stretching or di- 
agonal strain r\ in all three directions relative to the ref- 
erence substrate; i.e., u{(x) = i]Xi for i — x,y, z. 

This treatment assumes a fixed misfit at the interface, 
and would be exact for an interface between two crys- 
tals with a fixed epitaxial relationship, such as Si-NiSi2. 
However, the case of greatest interest, Si-SiC>2, is more 
complex. The large volume increase upon oxidation is 
largely accommodated by expansion normal to the inter- 
face, with only a modest residual misfit stress. Since the 
microscopic oxidation processes that determine the resid- 
ual stress are unknown, our assumption is reasonable but 
untested in the context of Si oxidation. 

To compute the stresses one must solve the elastic 
force-balance equation V- o= for both substrate and 
film, subject to the boundary conditions a — ► as 

z — > — oo (in the substrate), o -hs — (force balance 
at the upper surface), u s {xi) = u*{xj) (continuity of dis- 

placement at the interface), and a -hi —a -hi (force 
balance across the oxide-substrate interface) |I3J]. 

For the planar interface (a q = 0), the displacement 
and stress must vanish in the substrate. For the film, 
only the z-component of the displacement is nonzero: 
u{(z) = efZ, with e/ = 3K f r)/(\ f + 2p, f ). The only 

<->/ 

nonvanishing components of a are the diagonal ones 
parallel to the interface: o~l. x = a' yy = — 2yJ e/. 

For the planar film, the oxidant current must be uni- 
form throughout the oxide. So, by solving Eq. (^) with 
the boundary conditions (^), (0), one readily calculates 
the growth rate of the oxide thickness l(t) from Eq. ([[]): 

£]( f ) = ( c * - Pch/h)/pf (q] 
dt K ' l/k x + l/h + l{t)/D ' 1 ' 

This result is essentially identical to that of Deal and 
Grove |l0|, viz., l(t) grows as t at short times and as t 1 ! 2 
at long times. Note that in practice c* > k 2 pf /k\, so l{t) 
always increases with time. 

For the modulated geometry of Fig. 1, the interface 
position is Zi(x, t) — Zq{€) + a q (t)cos(qx), Zo(t) being its 
average position at time t. We calculate the curvature, 
displacements, and stresses to linear order in a q /\. With 
Eq. (^) one can then express SF t — 8Fi, and hence the 
oxidant flux j{xi) at the interface, in terms of a q , A, the 
oxide thickness I, and the parameters of the problem, 
such as 7, and erf'". One then solves V 2 c = in the ox- 
ide, to 0(a q /X), subject to the boundary condition (^) 
at the upper surface, and (||) at the interface. Knowing 
c(xi), one uses Eq. ([!]) to determine the interface posi- 
tion Zi(x,t): zo(t) is given by the solution of the planar 
problem, and the amplitude, ot q (t), of the perturbation 
evolves according to: 



da q /dt — Q q a q , (10) 

where 

Q q = Pf~ X (-w + wiq - w 2 q 2 ) ■ (11) 

The wo, wi, and W2 terms respectively represent the ef- 
fects of diffusion, stress, and interfacial tension. Since 
wo > and u>2 > 0, at sufficiently long or short wave- 
lengths fl q is always negative and perturbations decay 
with time. However, depending on the sign and magni- 
tude of the stress term wi , £l q may be positive for some 
range of q. Perturbations with wavelengths in this range 
will then grow exponentially with time, destabilizing the 
planar interface. 

The Wi can be written 

wo = -^(C Q ~ki-r c h)fo(ql) ; (12) 
D PS 

Wi = -f3e f l(Coh>Ji ~ r c k 2 o^)f^(ql) + 

(Coh^ - r c k 2 al )f*(ql)} ; (13) 
w 2 =P{C k l ll+r c k 2l2 )h{ql) . (14) 

Here 

_ c^pj 1 + r c k 2 /h + r c k 2 l/D 
° = (1 + k/h + hl/D) ' 

and r c = p c /pf is of order unity. The quantities fo(ql), 
fi(ql), ff(ql), and f 2 (ql) are complicated dimensionless 
functions of ql and the various parameters. The depen- 
dence of these functions on ql is weak and they are of 
order unity over the range of qVs of interest. For many 
reactions, including the oxidation of silicon, k 2 is negli- 
gibly small. Thus the reverse reaction (decomposition of 
the oxide) doesn't occur in practice. In this case, formu- 
lae simplify, since k 2 can be set to 0. 

The explicit g-dependence of the last two terms in Eq. 
( |lT| ) results from the strains (8e nn and 8e pp ) and the cur- 
vature (k) behaving like q and q 2 , respectively. Thus the 
smoothing influence of the oxidant diffusion term dom- 
inates and maintains stability at small q. This is the 
inverse of diffusion-limited growth, where the role of dif- 
fusion is destabilizing, as in the Mullins-Sekerka insta- 
bility. The interface tension term prevents perturbations 
of very small wavelength from growing, stabilizing the 
planar interface at large q. 

Thus planar growth is unstable when A = w\ — 
Awqw 2 > with the instability Q{q) > occurring only 
over the range of wavelengths g_ < q < q + . Here 



q± = {wi ± y wl - 4w w 2 )/2w 2 
This is illustrated in Fig. 2. 
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FIG. 2. An illustration of the stability calculation results. 
Notice that the interface is always stable (fl q < 0) in both the 
small and large q limit, and only become unstable (Q q > 0) 
in the intermediate range g_ < q < g+ when A > 0. 

It remains to determine the typical range of wavenum- 
bers q for which the instability holds. We first simplify 
the problem by assuming that ki is indeed negligibly 
small. As mentioned above, this assumption is typically 
well justified in practice. We must also estimate the pa- 
rameters a i and a\. Because oxidation is accompanied 
by expansion, the oxidation reaction should proceed more 
readily when the existing material near the interface has 
been dilated, i.e., when Se nn and 5e pp are positive. In 
other words, one expects both coefficients er" an d <j\ to be 
negative. Our elastic theory calculation shows, however, 
that f?(ql) < and f*(ql) > (with /f(<z0+/f (?Z) > 0), 
so the sign of w% is determined by the relative magnitudes 
of cr™ an d a i- Thus the instability presumably will not 
occur for all values of parameters. 

To see that it can nonetheless occur for some reason- 
able parameters, we estimate cr™ and a\ by assuming 
a i ~ a i ~ — 2/i/e', which makes W\ > 0. For the other 
parameters, we use rough numbers for wet oxidation in 
silicon jnj . Due to the large volume change during oxida- 
tion, the "effective" strain parameter e' is of order unity: 
e' ~ 1. Taking l(t) <C D /k\ (the most favorable limit 
for the occurrence of the instability), c* ~ 3xl0 19 cm -3 , 
k% ~ 5xl0 -5 cm /sec, D ~ 10~ 9 cm 2 /sec, pf ~ 10 22 cm -3 , 
W2 ~ 10 9 /cm — sec, and (3^ ~ 1.6xl0 24 cm -3 at about 
900°C, we find that the inequality W\ > 2^/wqw 2 , and 
hence the instability, holds for all e/ greater than about 
2 x 1CP 4 , which can indeed be achieved for modest 
values of misfit strain r\. For the numbers above, we 
find, roughly, g_ ~ 2xl0 4 /cm, and y+ ~ 5xl0 7 /cm for 
€f = 0.01. This implies an instability for wavelengths 
between roughly 1 nm and 0.3 /zm. 

Eq. (0) shows that the characteristic time tj over 
which the instability develops is set by l/^ 9max = 
<tw 2 pf/(wl - 4w wi), where q max = (q + + q~)/2 is the 
wavenumber for which a q grows most rapidly. For the 
sample parameters above, one finds that Tj *~ 1 sec. 
This time scale is short compared with typical process- 



ing times, or with the other characteristic time tx — 
Dpf/2c*k\ ~ 10 4 sec, which is roughly the time when 
the motion of the planar interface crosses over from t to 
i 1 / 2 behavior. 

Finally, we note that there is some intriguing exper- 
imental evidence of roughness in Si — SiC>2 interfaces 
for thin oxide layers |jj at 900°C. However, because of 
the poorly-understood complexities of real Si oxidation, 
with a volume expansion of order unity but a far smaller 
residual misfit, our results should be viewed as indicative 
of the kind of behavior that may occur, rather than as 
quantitatively applicable. 
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